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Most applications of time-dependent density-functional theory (TDDFT) use the adiabatic local- 
density approximation (ALDA) for the dynamical exchange-correlation potential V xc (r, t). An exact 
(i.e., nonadiabatic) extension of the ground-state LDA into the dynamical regime leads to a V^ c (r, t) 
with a memory, which causes the electron dynamics to become dissipative. To illustrate and explain 
this nonadiabatic behavior, this paper studies the dynamics of two interacting electrons on a two- 
dimensional quantum strip of finite size, comparing TDDFT within and beyond the ALDA with 
numerical solutions of the two-electron time-dependent Schrodinger equation. It is shown explicitly 
how dissipation arises through multiple particle-hole excitations, and how the nonadiabatic extension 
of the ALDA fails for finite systems, but becomes correct in the thermodynamic limit. 



I. INTRODUCTION 

The essential ingredient of time-dependent density- 
functional theory (TDDFT) HQ, the exchange-correla- 
tion (xc) potential V xc (r,t), is frequently obtained using 
the adiabatic local-density approximation (ALDA): 

v;t LDA (r,t) = %M , (1) 

where e xc (n) is the xc energy per particle of a homoge- 
neous electron gas of density n. The adiabatic approxi- 
mation means that all functional dependence of V xc (r,t) 
on prior time-dependent densities n(r',t'), t' < t, is 
ignored. Neglecting the retardation implies frequency- 
independent and real xc kernels in linear response. This 
approach has been widely used in quantum chemistry for 
calculating molecular excitation energies [3j- 

The adiabatic approximation is known to work well 
for excitation processes in many-body systems that have 
a direct counterpart in the Kohn-Sham system, such as 
atomic and molecular single-particle excitations. On the 
other hand, for more complicated processes such as dou- 
ble or charge-transfer excitations the ALDA can fail dra- 
matically 0, d| . A recent study [|| has shown that the 
ALDA can completely break down in dynamical pro- 
cesses where the electronic density rapidly undergoes 
large deformations. 

Several recent papers have addressed the question how 
the LDA for ground-state calculations should be p ro p erl y 
extended into the dynamical regime [EBB II EI El 01- 
Vignale and Kohn showed that a nonadiabatic lo- 
cal approximation for exchange and correlation requires 
the time-dependent current j{r,t) as basic variable (C- 
TDDFT). This formalism was later recast in the language 
of hydrodynamics, where xc effects beyond the ALDA ap- 
pear as viscoelastic stresses in the electron liquid [3, UJJ ■ 
An alternative nonadiabatic theory formulates TDDFT 
from the point of view of an observer in a co-moving La- 
grangian reference frame (L-TDDFT) [l^. In Ref. @, 
the technical details of C-TDDFT and L-TDDFT are 
critically examined and compared. 



To date, most applications of TDDFT beyond the 
adiabatic approximation take place in the frequency- 
dependent linear-response regime. A major success of C- 
TDDFT was the work by van Faassen et al. El wn0 cs ^~ 
culated static axial polarizabilities in molecular chains, 
with much improvement over the ALDA. 

The C-TDDFT formalism has recently been applied to 
describe linear and nonlinear charge-density oscillations 
in quantum wells by solving the time-dependent Kohn- 
Sham (TDKS) equation [14j. It was shown that the retar- 
dation caused by the memory of the xc potential has the 
striking consequence of introducing decoherence and en- 
ergy relaxation, i.e., the oscillating density experiences a 
damping. The mechanism causing this behavior has been 
discussed by D'Agosta and Vignale Technically, dis- 
sipation arises in C-TDDFT from a velocity-dependent 
xc (vector) potential which breaks the time-reversal in- 
variance of the TDKS Hamiltonian. As a result, a system 
tends to relax from a nonequilibrium initial state to an 
equilibrium final state with higher entropy. But where 
does the dissipated energy go? 

Because the system is closed and isolated (there is 
no coupling to a thermal bath), the total energy should 
be conserved. According to Ref. EI dissipation in C- 
TDDFT has to be understood in the sense that energy is 
redistributed between two subsystems with different sets 
of electronic degrees of freedom, coupled by Coulomb in- 
teractions. In the quantum well examples of Refs. El an d 
Il5l the transfer of energy occurs from a collective motion 
along the confinement direction into low-lying lateral ex- 
citations of the two-dimensional electron gas in the quan- 
tum well plane. However, in C-TDDFT this transfer pro- 
cess is never directly observed, since the TDKS equations 
are solved only for the electron dynamics perpendicular 
to the quantum well plane. 

The purpose of the present paper is to give an explicit, 
pedagogical illustration of the road towards dissipation 
in collective electronic motion. We will consider a two- 
electron model system that is simple enough so that its 
dynamics can be treated numerically exactly via solu- 
tion of the full time-dependent Schrodinger equation, and 
compare it with TDDFT within and beyond the ALDA. 



In particular, we will focus on charge-density oscillations 
along one direction of the system, and how the exact cal- 
culations show that the amplitude of these oscillations 
changes over time. This amplitude modulation comes 
from a superposition of transitions between the ground 
state and singly excited states and between singly and 
doubly excited states, including a coupling to the trans- 
verse degrees of freedom due to Coulomb interactions. 

In ALDA, all effects involving multiple excitations 
are completely absent; in C-TDDFT, multiple excita- 
tions are implicitly included, but for finite systems their 
contribution is strongly exaggerated, producing an un- 
physical damping. Based on the insights of our simple 
two-electron system, we will discuss how the dissipative 
behavior emerges in the thermodynamic (large-system) 
limit, and to what extent it is then correctly described 
by C-TDDFT. 

In section [H] we give the technical details of our 
two-electron model system and describe how the full 
Schrodinger equation and the TDKS equations with and 
without memory are solved. Section ITTTI gives our results 
and discusses the physical process of dissipation of col- 
lective charge-density oscillations. Conclusions are given 
in section ITVl 
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FIG. 1: (Color online) Schematic illustration of one cycle of 
a charge-density oscillation of a two-electron system on a 2D 
quantum strip of width A and length L. Darker areas rep- 
resent regions of charge accumulation. Snapshots are shown 
at times as indicated, where T is the duration of one cycle. 
The model assumes periodic boundaries along x and hard-wall 
boundaries at z = and z = A. 



II. MODEL SYSTEM 



Consider two electrons on a two-dimensional (2D) 
quantum strip of length L and width A, positioned in the 
x — z plane. In the following, we will be mostly interested 
in situations where L » A. The system has hard- wall 
boundary conditions at two ends of the strip, at z = 
and z = A, and periodic boundary conditions along the 
^-direction. In other words, the electrons are living on a 
strictly 2D surface whose topology is equivalent to that 
of a cylindrical tube of length A and circumference L. 

On this 2D quantum strip we first calculate the elec- 
tronic ground state in the presence of a linear external 
potential which depends only on z: 



V(z) = Fz , 



(2) 



where F is a constant field strength. At the initial time 
t = 0, this external potential is suddenly switched off, 
which triggers a charge-density oscillation along z (see 
Fig. 0. The electronic density thus remains uniform 
along the x-direction for all times. The goal is to fol- 
low the time evolution of the system for many cycles 
of the charge-density oscillations, comparing the exact 
numerical solution of the two-electron Schrodinger equa- 
tion with TDDFT solutions within and beyond ALDA. 
Atomic (Hartree) units are used throughout. 



A. Two-electron Schrodinger equation 

1, Ground state 

The static two-electron Schrodinger equation for our 
problem reads 
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-T-T + ^ + ^ + , ri -r 2 | 
x$,-(r 1 s 1 ,r 2 s 2 ) i 



(3) 



where ri )2 = (xi^, 21,2), and si j2 denotes the spin. We 
expand the two-electron eigenstates $j in a basis of Slater 
determinants: 



^•(risi,r 2 s 2 ) = C U 



where 

^= fi I/o K, 



(4) 
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ip„ 2 k 2 [xi, zi)£2(si)VWfci ( x 2, z 2 )Ci(s 2 ) ■ (5) 



Here, £ are single-particle spinors, and for the spatial 
part we choose the non-interacting single-particle wave 
functions for constant external potential: 



ipnk(x, z) 



LA 



e lkx sinnz 



(6) 
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with quantum numbers 
2ttk 



K = 0,±l,±2,. 



n= — , ^ = 1, 2, 3, . 



(7) 

(8) 



In other words, we satisfy the given boundary conditions 
(see Fig. ^) by choosing plane wave basis states along 
the strip and standing waves across the strip. 

Inserting the basis expansion (@J into the Schrodinger 
equation yields the following equation for the expan- 
sion coefficients: 



2. Time evolution 

Once equation © has been diagonalized, the next step 
is to determine the time evolution of the ground state 
$i(r!Si, r 2 s 2 ,t) after the linear external potential has 
been switched off. Rather than explicitly solving the 
time-dependent two-electron Schrodinger equation, this 
is most easily done by expanding $i in the complete set 
of field-free eigenstates, defined as follows: 



E 



${ = 0. 
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Thus, 



^V\V 2 K\K 2 J 



(9) 



* x (t)=x;^(*)*/ 



(15) 



(16) 



(17) 



Here, the kinetic energy and external potential matrix 
elements are given by 



where 



T — 



n 2 + k 2 



(10) 



A 3 (t) = expl-iEjt] Ci'/ V2KlK2 Cl ll/2KlK2 . (18) 
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and 



From this, we obtain the time-dependent density as fol- 
lows: 
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Vfiv = -r / dz sinmzsinnz V(z) . (11) r 

AJo n(z,t) = / d27 -2|*i(risi,r 2 s 2 ,i)| 



The matrix elements of the Coulomb interaction are 

4 f A f A 

Wgfigg = dz x dz 2 sin To! Z! sinTO 2 z 2 

x sinnizi sin n 2 z 2 S kl +k 2 , qi +q 2 Ik 2 -q 2 {zi, z 2 ) 

(12) 
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^1^2 
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+ sinm 2 z sinn 2 z ^ ^] , 



(19) 



where 



where 



I k -q(zi,Z 2 ) 



dx 



cos[(fc — q)x] 



^/x 2 + (zi - z 2 ) 2 

= 22f [|fc-g||2!i-«2|] ) M<7 (13) 
= -21og|z 1 -z 2 |, fc = g. (14) 

Here, Kq is a complete Bessel function of the second kind 
in standard notation, and in the case k = q an additional 
divergent term is cancelled by the positive background. 

Solving equation jnj numerically one finds that a rela- 
tively small basis size including states with with no more 
than k = ±10 and v = 10 [Eqs. Q, JSJl] is sufficient. 
The computational task is therefore quite manageable. 

Furthermore, it turns out that, due to symmetry and 
momentum conservation, only those two-electron basis 
states ^u 1 v 2 k 1 k 2 contribute which have zero net current 
along the strip, i.e., only states with K4 = — n 2 are 
needed. This corresponds to two-electron states where 
one electron travels to the right and the other to the left. 



Finally, the time-dependent dipole moment is 



d(t) = dz zn{z 1 1) 
Jo 



(20) 



(21) 



B. TDDFT 

1. Ground state 

The two-electron problem described above can be 
solved, in principle exactly, using the TDKS formalism. 
We begin with the static Kohn-Sham (KS) equation: 



-— + V(z) + V H {z) + V xc (z) - E n 



b n {x, z) = . 

(22) 
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This equation separates in x and z, and we make the 
ansatz 



4>n{x, z) 



E n — 



1 

— 7= e 

h 2 k 2 

2m 



-3 i 



(23) 
(24) 



where the index k is given by Eq. (JJJ). The ground-state 
solution has k = 0, and we end up having to solve a 
one-dimensional equation for <fj(z) and £f. 



2dz 



2 +V{z) + V a (z) + V„{z)-e j 



(25) 

To solve the single- particle KS equation, we expand in a 
standing- wave basis as follows: 



N 



Cl sinnz , 



(26) 



where 



E 



2 "flu T K/iK T V^j/ T l 7 ^,, 



= . (27) 



From this, the ground-state density follows as 
4 



n(z) 



LA 



C X *C\, sin mz sin? 



(28) 



//// 



The matrix elements for the external, Hartree and xc 
potential are calculated from Eq. JUJ. The Hartree 
potential is given by 



V H (z) = -2 f dz'n(z') \og\ 



(29) 



plus a diverging constant which is cancelled by the 
positive background. For the xc potential we use 
the LDA within the parametrization of the 2D elec- 
tron gas of Tanatar and Ceperley 0. For the (spin- 
unpolarized) systems under consideration, the more mod- 
ern parametrization by Attaccalite et al. |17| gives al- 
most identical results. 

Figure [21 shows the two-electron ground-state density 
n(z) on a 2D quantum strip of width A = 10 and length 
L = 50, in the presence of a linear external potential @ 
with field strength F = 0.02. For these system param- 
eters, the 2D Wigner-Seitz radius r s = (7m)" 1 / 2 has a 
value of r s = 6 at the maximum of the density distri- 
bution. The agreement between the exact and the LDA 
density is reasonably good, and in fact becomes better 
for smaller quantum strips where the density is higher. 
In general, the LDA system is found to be a little more 
polarizable than the exact system. 



0.01 



0.008 




FIG. 2: Two-electron ground-state density n(z) on a quan- 
tum strip of width A = 10 and length L — 50, with field 
strength F = 0.02. Full line: exact solution. Dashed line: 
LDA. 



2. Time evolution 

As for the two-electron case, we set the charge-density 
oscillations in motion by suddenly switching off the ex- 
ternal potential at the initial time t — 0. The task is to 
solve the TDKS equation 



d_ 

Of 



<p{z,t) = 0, (30) 



with initial condition <p(z,0) = tp\(z). The time- 
dependent KS orbital tp(z,t) is expanded similar to 
Eq. H26[l. an d the time-dependent expansion coeffi- 
cients Cfi (t) are numerically determined using the Crank- 
Nicholson algorithm plus predictor-corrector scheme [2j. 

In the following, we will consider V xc (z,t) within and 
beyond the ALDA. The C-TDDFT expression for a time- 
dependent xc potential with memory is called ALDA+M, 
and written as follows 0, 0] : 



V xc (z,t) = V£ WA (z,t) + V™(z,t), 



(31) 



with the memory part 



dz' 



V z 'f7 XC:ZZ /(z', t) 



l n(z',t) 

The zz component of the xc stress tensor is given by 



(32) 



,(z',t) 



Y(n(z', t),t- t')\7 z ,v z ,{z', t)dt' . (33) 



Here, v(z,t) = j(z,t)/n(z,t) is the time-dependent ve- 
locity field, where j(z,t) is the current density. In 2D, 
the memory kernel Y is given by 

r(M-t')=fec-- / — 3/ x L c Mcos[o^-0], (34) 

7T J UJ 

with the 2D xc shear modulus of the electron liquid 

Mxc = n 2 (9f?/ x L c (0) - (ne xc )") (35) 
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(the prime denotes a derivative with respect to n). In 
the following, we use the Holas-Singwi parametrization 
for the longitudinal frequency-dependent xc kernel of the 
2D electron liquid [l9| : 



Alj 



B 2 



AB 



B 2 +LU 2 

The coefficients A and B are given by 

IItt 2 



A = — 



32 



B 



A 



with 



/oo — 



2n 



(ne xc )" - f a 




(36) 
(37) 

(38) 
(39) 

(40) 



It is easy to see that this simple parametrization for f£ c 
leads to zero shear modulus, /i xc = 0. A more sophis- 
ticated interpolation formula, with finite /i xc , has been 
derived by Qian and Vignale 01 > but its input param- 
eters are currently only available for a limited range of 
densities in the metallic regime. For our purposes, it is 
therefore preferable to work with the Holas-Singwi for- 
mula Ij36(l . which has the additional advantage that it 
leads to a very simple expression for the memory kernel: 



Y(n,t-t') 



An 2 



-B(t-t') 



(41) 



i.e., the system experiences an exponential memory loss. 
This is similar to what was observed [6J, Il4j in 3D sys - 
tems using the Gross-Kohn parametrization for f^ c |20| . 
As a consequence, the numerical evaluation of the time 
integral in Eq. I|33|) can be simplified by introducing a 
cutoff in t — t', i.e., not the entire history of the system 
from t — onwards needs to be included. 



III. RESULTS AND DISCUSSION 
A. Charge-density oscillations 

Figs. |3] and 2] compare the time-dependent dipole mo- 
ment d(t) [Eq. pijl ] for quantum strips of width A = 10 
and lengths L = 50 and L — 100, respectively, calculated 
from the exact density l|19(l and from the ALDA density. 
The initial state was prepared with an external poten- 
tial @ of field strength F = 0.02, which was abruptly 
switched off at t = 0. 

At first sight, the ALDA charge-density oscillations 
seem to agree well with the exact ones, as far as the fre- 
quency and the average amplitude of d(t) are concerned. 
On closer examination, however, we observe a beating 




4000 



FIG. 3: Time-dependent dipole moment d(t) associated with 
the charge-density oscillations in a quantum strip of width 
A = 10, length L = 50, and initial field strength F = 0.02. 
Top: exact solution. Bottom: ALDA. 
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FIG. 4: Same as Fig. but for a strip of length L = 100 



pattern in the exact charge-density oscillations, which 
shows up as a low-frequency modulation of the ampli- 
tude of d(t). This effect is not reproduced by the ALDA, 
which produces a constant amplitude (we disregard here 
the small, rapid wiggles in the amplitudes of d(t), which 
are a nonlinear effect caused by the relatively strong ini- 
tial field, see Fig. In the following, we will focus on 
discussing the origin of these modulations, and on the 
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TABLE I: Leading terms in the expansion 117H of the time- 
dependent two-electron wave function $(£) in terms of field- 
free eigenstates (for a strip with A = 10, L — 100, and 
F — 0.02). (ui,vs) indicates the dominating single-particle 
configurations (see Fig. |SJ. 



J 




A) 


Ei 


1 


(1,1) 


0.897406 





2 


(1,2), (2,1) 


0.098766 


0.148661 


3 


(2,2) 


0.002698 


0.294534 


4 


(1,3), (3,1) 


0.001008 


0.394043 


5 


(2,3), (3,2) 


0.000051 


0.542683 


6 


(1,4), (4,1) 


0.000048 


0.739263 



related shortcomings of the ALDA and the consequences 
thereof. It will turn out that this effect provides a key 
to understanding the meaning of dissipation in TDKS 
theory. 

In section lll A 2l we considered the time evolution of the 
exact two-electron state, and explained how it can be ob- 
tained by expanding the time-dependent wave function 
in the complete set of field-free eigenstates, see Eq. ((T7I) . 
For the system parameters and field strengths under con- 
sideration, this expansion turns out to be dominated by 
just a few leading terms in the summation over Aj(t)$j. 
Looking at those few terms with the largest |Aj(i)| 2 will resentation $ 



give us sufficient information to understand the electron 

dynamics leading to the beating pattern in d{t). 

Let us analyze in detail the case A = 10, L = 100 and 

F = 0.02. We have solved the two-electron Schrodinger 

equation with 11 plane- wave and 8 standing- wave basis 

states, i.e., \k\ < 5 and v < 8 in Eqs. and ©. Tabled 

shows the six leading terms in the expansion l|17|) of $(i) 

in terms of field- free eigenstates, i.e., those terms with 

the largest Aj , and the associated energies Ej , where we 

define the field- free ground-state energy to be E± = 0. 

The remaining terms in Eq. i|17fl have values of A| that 

are orders of magnitude smaller. 

According to equation J2J, each field- free eigenstate 
f 

$j is represented as a sum of single-particle Slater deter- 
minants. The leading <I>J's are dominated by configura- 
tions ^u 1 v 2 k 1 k 2 whose standing-wave quantum numbers 
{v\,V2) are given in the second column of Table [I] and 
which can have a broad range of plane waves (ki,K2)- 
This is illustrated in detail in Fig. \5\ which shows his- 

in the expansion 
of the first four lead- 



tograms of the coefficients [C; 



/ _ 



V\V 2 K\K\ 



S "1 

ing field-free eigenstates. One can clearly see that there 
are dominating pairs of standing- wave quantum numbers 
iyit v 2), which explains the assignment in the second col- 
umn of Table [I] Each configuration with standing waves 
V2) along z is accompanied by left and right running 
waves (k, —k) along x. The case n — is dominant, but 
finite re are not negligible. 

To understand the beating pattern in the dipole oscil- 




FIG. 5: Values of the coefficients \Ci^ V2K1K ^ 2 in the rep- 



of the first four 



^ i^l C[/^ 2 kiK2 "l^^l-: 

leading field-free eigenstates (see Table[IJ. Each configuration 
with standing waves (fi,i*j) along z also has running waves 
along x with (re, —re), where —5 < re < 5, shown here as his- 
tograms. The dominating configurations are for re = 0, but 
finite re are not negligible. 



lations of Figs. OJand^l we now focus on the first three 
leading field-free eigenstates and their standing-wave 
quantum numbers (1*1, and for the moment disregard 
the running waves along x. The beating pattern essen- 
tially arises from a superposition of two dipole oscillations 
associated with the transitions (1, 1) — > (1, 2), (2, 1) and 
(1,2), (2,1) — > (2,2). The associated energy differences 
are E 2 — E± = w 2 i = 0.148661 and E 3 — E 2 = uj 32 = 
0.145873. The two oscillation frequencies u>2t and W32 
are very close, and their difference oj 2 i — <^z 2 — 0.002788 
is precisely the frequency of the amplitude modulation 
of d(t). The resulting modulation period is T moc i = 
27r/(w2i — W32) = 2254, which agrees extremely well with 
the data shown in Fig. 0] Similarly, for the case L = 50 
shown in Fig. [31 we find T mo( j = 964 (here, the difference 
W21 — w 32 is a bit bigger). The amplitude of the modula- 
tions of d(t) depends on the field strength F and remains 
small as long as A\,A\ ^> A 3 . 

It is now easy to see why the ALDA misses the beat- 
ing pattern in d(t): the reason is that it does not ac- 
count for doubly-excited configurations. The ALDA in- 
cludes only single excitations, which are the only possible 
excitations of the KS system. Thus, transitions involv- 
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ing the (^1,^2) = (2,2) configuration, which are crucial 
to explaining the modulation of d(t), cannot occur, and 
therefore no superposition effect takes place. 

In addition to the standing- wave double excitations, 
the contribution of doubly excited running-wave states 
(k, —k) along x are also important. Again, the ALDA 
only includes the case n = (single excitations along x 
are not possible due to momentum conservation). On 
the other hand, ignoring the states (k, —k) with finite k 
in the expansion of the full two-electron wave function 
would lead to substantially different energies Ej , and the 
low- frequency beating pattern of d(t) would be destroyed. 

The exact xc potential (which we will construct in the 
next subsection) has to compensate for the absence of 
multiple excitations in the TDKS wavefunction, and it 
does so through a nonadiabatic contribution. This is 
known from linear response theory Q , where the xc ker- 
nel must have a frequency dependence to describe double 
excitations. 



B. Exact xc potential and time-dependent energy 

1. Construction of the exact time- dependent xc potential 

If the density n(r,t) of a system of two electrons in a 
singlet state is given, it is a straightforward affair to con- 
struct that xc potential V xc (r,t) which, when employed 
in a TDKS equation, reproduces this density |2l|. The 
doubly occupied TDKS orbital can be written as 



(42) 



where the phase a is a real function and related to the 
current density as follows: 



Va(r,t) =j(r,*)/n(r,i) . 



(43) 



Inserting the ansatz l|42|l into the TDKS equation, one 
obtains 

V xc (r,t) = (r,t) + K d c yn (r, t) . (44) 
The first term, 

K S c at M) - iv 2 lnn(r,t) + i|Vlnn(r,t)| 2 



V(T,t)-Vu(v,t) , 



(45) 



is identical to the expression for constructing the static xc 
potential from a given static two-electron density (2^ , ex- 
cept that all quantities are now taken as time-dependent. 
The second term, 

^ d / n (r,t) = -d(r,t)-l|Va(r,t)| 2 , (46) 

has no static counterpart and is therefore a truly dynam- 
ical contribution. 



2. The time- dependent energy 

Let us now consider the time dependence of the to- 
tal energy E(t) of the two-electron system. Since we 
are dealing with free charge-density oscillations of a fi- 
nite 2D quantum strip, i.e., there is no time-dependent 
external force, the total energy of the full many-body 
system must obviously be constant. This is easy to see 
for the exact time-dependent two-electron wave function 
$(f): according to equations (|T7Jl and JTSJl, we simply 
have E(t) = Y^j Ej for all times. 

On the other hand, it is less obvious what the TDDFT 
expression for the exact time-dependent total energy of 
a many-body system should be (although some compo- 
nents of the TDDFT energy have been studied in Ref. 
[2^1. However, for the purpose of this paper it is suf- 
ficient to define a quantity which we call the adiabatic 
energy, E a (t). For a two-electron KS system with a dou- 
bly occupied single-particle orbital, we have |15| 



E a (t) 



2 J dr <p*(r,t) 



1 A 

2 7 + A * 



V(r,t) 



<p{T,t) +Ea[n{t)] +E xc [n{t)} -(47) 



Here, En[n] is the Hartree energy functional, E KC [n} is the 
ground-state xc energy functional that was used in the 
calculation of the initial state for the TDKS time propa- 
gation, and A xc is the nonadiabatic xc vector potential. 
Clearly, E a (t) is not the true energy of the many-body 
system, but it reduces to the ground-state energy in the 
static limit. It can then be shown that the rate of change 
of the adiabatic energy is [l5l| 



K(t) 



dr j(r,t) • A xc (r,i) 



(48) 



In the case of our 2D quantum strip, the density is spa- 
tially inhomogeneous along the z direction only, which 
means that we can replace the xc vector potential by the 
dynamical scalar potential, using the relation 



A xc (r,i) = -W x d c yn M). 



(49) 



Furthermore, the z-components of the physical current 
j z {z,t) and the KS current jKS,z(z,t) = 2 < 3(ip* dip / dz) 
of this particular two-electron system become identical. 
Therefore, we obtain 



Kit) 



-L I dz Jz (z,t)±V x d r(z,t). (50) 







dz 



This shows that the rate of change of the adiabatic energy 
is determined by the work done by the forces associated 
with the dynamical xc potential. 

In the case where we start from the exact time- 
dependent density and ask what the associated exact adi- 
abatic energy is, a direct evaluation of expression (|47|l is 
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FIG. 6: Full lines: cycle- averaged adiabatic energy E a (i) for 
the exact solutions of the charge-density oscillations on the 
2D quantum strip with L — 50 (top) and L = 100 (bottom). 
Short and long dashed lines: E a (t) calculated with ALDA 
and ALDA+M. For clarity, all energies are shifted so that 
they initially coincide. 



produces any (reasonably well-behaved) time-dependent 
density n(r,t) of a many-body system. In particular, if 
the exact functional for the xc potential V xc (r, t) is used, 
the TDKS system gives the exact time-dependent den- 
sity. However, the adiabatic energy (as defined above) is 
not required to remain constant like the true energy of 
the many-particle system (in the absence of an external 
field, of course). 

How does this happen in our 2D model system? The 
exact TDKS system has to somehow reproduce the am- 
plitude modulations of the time-dependent dipole mo- 
ment, but it cannot do so through a simple superposition 
of oscillations associated with single and double excita- 
tions, as it happens in the full two-electron Schrodinger 
equation: there are no double excitations in the TDKS 
system. Instead, the TDKS system has to produce the 
beating pattern in d(t) through the action of the xc po- 
tential. In other words, the dynamical part of the xc po- 
tential, V^p 1 , acts in a sense like an "external" potential 
which periodically drives and damps the charge-density 
oscillations of the system in order to increase or diminish 
the amplitude of d(t) . From this point of view, V^ yn thus 
periodically dissipates energy from the KS system, and 
then pumps it back into it. The adiabatic energy E a (t) 
serves as an indicator for this behavior. 



C. Memory effects and dissipation 



not possible, since one doesn't know the form of the exact 
xc energy functional. Fortunately, the adiabatic energy 
can easily be obtained from Eq. I|5U|I through a sim- 
ple time integration: E a {t) — J Q E a (t')dt' . This is a very 
convenient way of determining the exact adiabatic energy 
(up to an irrelevant constant) from the exact density. 

We have calculated E a (t) for the exact solutions of 
the charge-density oscillations for the 2D quantum strips 
with L = 50 and L = 100, see top parts of Figs. |3]andE] 
As expected, the exact adiabatic energy is not constant, 
but rather rapidly fluctuating with time. For the sake of 
clarity, and since these rapid fluctuations are not what 
we are primarily interested in, we define a cycle-averaged 
adiabatic energy E a (t), which is obtained by averaging 
the adiabatic energy at each time t over a time window 
of one period of the charge-density oscillation (duration 
- 40 a.u.). 

Figure shows the adiabatic energy for the two quan- 
tum strips, calculated with different methods [both the 
ALDA and ALDA+M energies have been directly de- 
termined from Eq. (|17Jl : the ALDA+M results will be 
discussed in the next subsection]. One clearly sees slow 
oscillations of the exact E a (t) with the same period as the 
amplitude modulations of d(t) (see Figs. |3]and0J|. By 
contrast, E a (t) in ALDA is constant as expected, since 
the ALDA dipole amplitudes are not modulated. 

These results provide some interesting insights into 
TDDFT. It follows from the Runge-Gross theorem [| 
that there exists a unique TDKS system which re- 



in the previous subsections, we have seen that the 
ALDA fails to reproduce some key features of the dynam- 
ics of the two-electron system, related to the fact that it 
misses the double excitations. We have also seen that 
the exact xc potential makes up for this through nona- 
diabatic contributions. Let us now see how the nonadia- 
batic ALDA+M approximation of C-TDDFT performs, 
which was described in Section III B 21 

Figure \7\ shows the time-dependent dipole moment 
for the charge-density oscillations of the same 2D quan- 
tum strips discussed above, but now calculated using the 
nonadiabatic ALDA+M xc potential [Eq. (|5T|) ] . In both 
cases, the dipole oscillations are exponentially damped, 
similar to what was previously observed for plasmon os- 
cillations in a semiconductor quantum well [14|. Repre- 
senting the dipole moment as d(t) ~ do cos(ujt)e~ rt , we 
find a damping rate of about T = 0.0007 for the L = 50 
strip and T = 0.00055 for the L = 100 strip. As ex- 
pected, the system with L — 100 has a somewhat weaker 
damping, due to its lower particle density, which reduces 
the probability of electron-electron scattering. 

The long-dashed lines in Fig. 0show the adiabatic en- 
ergy E a (t), calculated with ALDA+M. In both quantum 
strip systems, the energy is dissipated at an exponential 
rate, E a (t) — E a (0)e~ 2Tt . As discussed in the previous 
subsection, the exact adiabiatic energy oscillates, but on 
average there is no dissipation. These results clearly show 
that the nonadiabatic ALDA+M functional fails for our 
finite two-electron system. It results in an unphysical 
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FIG. 7: Time-dependent dipole moment calculated with 
ALDA+M, for the same charge-density oscillations treated 
in Figs. [3]and0] Upper panel: L — 50, lower panel: L = 100. 

damping of the charge-density oscillations. 

Technically, the dissipation arises from the fact that 
V™ is a velocity-dependent potential. As was discussed 
in Ref. Q, the history-dependence of V™, which is gov- 
erned by the memory kernel Y(n,t — t'), is such that it 
accounts for both dissipative and elastic properties of the 
electron liquid. The microscopic properties of the many- 
body system enter through the frequency-dependent xc 
kernel of the homogeneous electron gas, / xc (w). This 
function describes dynamical processes of the homoge- 
neous electron gas that go beyond single-particle ex- 
citations, i.e., / xc (w) contains the physics of multiple 
particle-hole excitations. In principle, this should be a 
good thing, since we have seen from our model that the 
main defect of the ALDA is the absence of double excita- 
tions. But in spite of this, the ALDA+M does not work 
in our system. How can one understand this? 



D. Thermodynamic limit 

The central reason for the failure of the ALDA+M for 
finite systems is the fact that it is based on the homoge- 
neous electron g cLS, 1.6.. cL reference system of infinite ex- 
tent. Since it is a local functional, the potential V^(z, t) 
has exactly the same value at some position x of our 
quantum strip as it would have if the strip were infinitely 
long, but with the same electron density n(z,t). Thus, 
by construction, ALDA+M treats all systems locally as 
if they were infinite, even if they're not. 

We can carry out the thought experiment of increasing 
the length L of the strip, and simultaneously adding elec- 



trons to keep the same n(z, t) as for the two-electron sys- 
tem. The ALDA and the ALDA+M xc potentials would 
be unchanged, likewise the Hartree potential, and we 
would find exactly the same charge-density oscillations 
across the strip. In particular, the damping in ALDA+M 
would be the same, irrespective of the length of the strip. 

On the other hand, the dynamics of the exact many- 
body system will change dramatically if we increase both 
L and the particle number. The time-dependent many- 
body wave function will contain not only single and dou- 
ble excitations, but a vast number of multiple excita- 
tions. The density of levels in the excitation spectrum 
will grow, and eventually turn into a continuum. Re- 
call that for the two-electron system, we explained the 
periodic amplitude modulation of the charge-density os- 
cillation through a superposition of two frequencies as- 
sociated with the dominant single and double excitation. 
If the system size grows, many more such transitions will 
play a role, and we will have to form a coherent super- 
position of many close-lying oscillators. The resulting 
beating pattern will become more complex, and seems 
difficult to predict. 

However, we can get a clue from comparing the case of 
L = 50 and L — 100. For the longer strip, the modulation 
period increases, i.e., the recurrence time becomes longer. 
In the limit of infinite system size, this suggests that the 
recurrence time will in fact become infinite. In other 
words, the charge-density oscillations will be irreversibly 
damped. 

But where does the energy go? In the free charge- 
density oscillations considered here, the wave function 
can be expressed as a linear superposition of field- 
free many-body states, each of which carrying a time- 
dependent phase exp(-iEjt) [this is a generalization of 
Eq. JT7J) to N particles]. Thus, the energy is, from the 
very beginning when the oscillation is triggered, shared 
in a fixed manner among all excited-state configurations 
that make up the time-dependent many-body wave func- 
tion. In turn, the charge-density oscillation is a coher- 
ent superposition of a continuum of single and multiple 
particle-hole excitations, which steadily run out of phase. 
This reduces the amplitude of the collective mode due to 
destructive interference. 

To give a simple illustration, consider the case of an ex- 
ponential damping, d(t) = d cos(wi)e _r *. We can carry 
out a Fourier analysis of the spectral content of d(t), and 
the result is that d(t) arises from a superposition of a con- 
tinuum of oscillators whose frequency distribution has a 
Lorentzian shape of half- width T centered around uj. 

The ALDA+M functional of C-TDDFT, per construc- 
tion, becomes exact in the limit of an extended system 
whose ground-state density, as well as the inhomogeneity 
of the time-de pen dent perturbation, are slowly varying 
in space 0, 1^. llOj. For such a system, the damping of 
the plasmon amplitude will be correctly described. In 
the many-body system, the damping occurs through in- 
terference within the continuum of multiple excitations, 
but in C-TDDFT, it has to happen in a completely dif- 
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ferent way: the nonadiabatic piece of the xc potential 
acts like an external damping force. The outcome, i.e. 
the behavior of the time-dependent density n(z, t), is the 
same. 



IV. CONCLUSION 

We have presented a simple two-electron system which 
has the appealing feature of being exactly solvable with 
modest computational cost, and rich enough to provide 
new insight into the problem of dissipation of collective 
electron dynamics. From the point of view of the exact 
time-dependent many-body wave function, plasmon dis- 
sipation occurs through a superposition of a continuum of 
oscillators, associated with transitions between multiply 
excited states, which slowly and irreversibly run out of 
phase. The phenomenon can be viewed like a beat, but 
with an infinitely long recurrence time. Consequently, 
there is no loss of energy in the many-body system, and, 
in a sense, not even a re-distribution into other degrees 
of freedom. 

From a TDDFT point of view, all we can say is that 
we have a time-dependent density which produces a time- 
dependent dipole moment whose amplitude steadily de- 
creases. The exact TDKS system accomplishes this 
through a nonadiabatic xc potential which acts like a 
damping force. As a result, energy of the KS system is 



lost, but this is the price we have to pay to reproduce the 
exact density. 

The ALDA+M xc functional of C-TDDFT has been 
constructed for infinite systems, and becomes exact in 
the appropriate limits For finite systems, on the other 
hand, it introduces a spurious damping of electron dy- 
namics. For example, if the method is applied to atoms, 
one obtains excitation energies with finite linewidths |24| . 
On the other hand, the static limit of the Vignale-Kohn 
functional B seems to work well for polarizabilities of 
polymers [I^- Thus, more tests of C-TDDFT are needed 
to explore its usefulness for practical applications. 

However, it seems unlikely that a time-dependent xc 
functional based on the homogeneous electron gas can 
correctly describe the subtle aspects of the dynamics 
of both finite and extended systems that we have dis- 
cussed in this paper. A more promising approach may 
be through orbital-based functionals such as the time- 
dependent optimized effective potential p5|. which will 
be the subject of future studies. 
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